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Abstract. Fortran subroutines to calculate helicity amplitudes with goldstinos, which appear as the longi- 
tudinal modes of massive gravitinos in high energy processes, are added to the HELAS (HELicity Amplitude 
Subroutines) library. They are coded in such a way that arbitrary amplitudes with external goldstinos 
can be generated automatically by MadGraph, after slight modifications. All the codes have been tested 
carefully by making use of the goldstino equivalence theorem and the gauge invariance of the helicity 
amplitudes. Hadronic total cross sections for associated gravitino productions with a gluino and a squark 
are also presented. 



1 Introduction 

Goldstinos are Goldstone fermions, massless spin- 1/2 par- 
ticles, associated with spontaneous supersymmetry (SUSY) 
breaking, and appear as the helicity ±1/2 states of massive 
gravitinos via the super-Higgs mechanism in local super- 
symmetric extensions to the Standard Model (SM). While 
the interactions of the helicity ±3/2 components of the 
gravitino are suppressed by the Planck scale, those of the 
helicity ±1/2 components are suppressed by the SUSY 
breaking scale and can be important even for collider phe- 
nomenology in low-scale SUSY breaking scenarios, e.g., 
gauge- mediated SUSY breaking [T]. 

In the recent paper [2] K. Hagiwara and the authors 
introduced new HELAS subroutines [J to calculate helic- 
ity amplitudes with massive spin- 3/ 2 gravitinos. They are 
coded in such a way that arbitrary amplitudes with exter- 
nal gravitinos can be generated automatically by MadGraph 

01- . 

In this paper, taking into account high energy pro- 
cesses with the center-of-mass (CM) energy -y/i ^ ^13/2, 
we present new HELAS subroutines for goldstino interac- 
tions based on the effective Lagrangian below, and imple- 
ment them into MadGraph/MadEvent(MG/ME)v4 [1115115], as 
an alternative to the code for gravitinos |2]Q As we will 
see later, in the high energy limit, the new code for gold- 
stinos agrees with the code for gravitinos with the correc- 
tion of order 7743/2 /%/« due to the goldstino equivalence 
theorem. The new code could be useful especially for col- 
lider phenomenology, where the goldstino limit is good 
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approximation for most of the cases. Practically, calcula- 
tions of helicity-summed amplitude squared for goldstino 
processes are faster than those for gravitino roughly by 
a factor of four due to the number of the helicity states 
and the simpler structures of the HELAS subroutines. We 
also note that the goldstino code can be applied to models 
such as broken global SUSY and goldstini Uj. 

The effective interaction Lagrangian for a goldstino in 
non-derivative form is [8u9ul0j 



PI ™3/2 



4V6A/pi TO3/2 

+ Jt"^' ^75A^")°0'*r(")>^- (1) 

V6A'/pim3/2 

with the reduced Planck mass il/pi = A/pi/a/Stt ^ 2.4 x 
10^* GeV and the gravitino mass 7713/2- is the Majorana- 
spinor goldstino field, and 0* are spinor and scalar 
fields in the same chiral supermultiplet, and Pr/l = ^(li 
75) are the chiral-projection operators. y("=3,2,i)a ^^iq 
SU{3)c (a = !,••• ,8), SU{2)l (a = 1,2,3) and U{1)y 
generators, respectively, and (?a=3,2,i are the correspond- 
ing gauge couplings. The field-strength tensors for each 
gauge group are 



p(3)a — a Aa _ 



d„B„, 
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and the corresponding gauginos A*^"^'^'^^^'" are gluinos 
winos {W"") and bino (B), respectively. 
The equivalence of the non-derivative form in ([Ij and 
the derivative form which is obtained by the replacement 
of the gravitino field ~ ^2/39^V/^3/2 in the grav- 
itino interaction Lagrangian (see, e.g. eq. (2) in [2]) has 
been proved in [5]. The following features of the interac- 
tion Lagrangian ([T]) are worth noting: (i) The ip-f-cjy-A^ 
vertex is absent, while a new quartic vertex, -f/^-A-i^i-^, ex- 
ists, (ii) The couplings are proportional to the mass split- 
ting inside the supermultiplet, m^; — rriji and m\, and 
inversely proportional to the SUSY-breaking vacuum ex- 
pectation value through the gravitino mass 

^3/2 = {F)/VSMpi. (5) 

The paper is organized as follows: In Sect. [2] we test 
our code by making use of the goldstino equivalence the- 
orem, and in Sect. [3] we give hadronic total cross sections 
of associated gravitino productions with a gluino and a 
squark as sample numerical results. Sect. 2] presents our 
brief summary. In App. |2]we give the new HELAS subrou- 
tines for goldstinos, and in App. |B] we describe how to 
implement the amplitudes into MG. 

2 Checking the goldstino-gravitino 
equivalence 

In this section, we check the new HELAS subroutines and 
the modified MG for goldstinos by using the goldstino equiv- 
alence theorem. As mentioned in introduction, goldstinos 
appear as the longitudinal modes of massive gravitinos 
and their interactions become dominant over the trans- 
verse modes in high-energy processes. By using the MG/ME 
package with gravitinos [2] , we test an agreement between 
the goldstino amplitudes and the gravitino amplitudes in 
the high-energy limit. 

We consider the goldstino production processes asso- 
ciated with a squark 

q + g q + G and q + g q + G, (6) 

which involves the 4'~f~4' vertex as well as ^-X-A^, and 
with a gluino for the gg initial state 

g + g^g + G, (7) 

which involves the ^p-X-Af^ and tp-X-Af^-A^, vertices. The 
Feynman diagrams shown in Fig. [T] and the correspond- 
ing helicity amplitudes are generated automatically by the 
modified MG. The details of the HELAS subroutines for gold- 
stinos and those implementation to MG are presented in ap- 
pendices. We note that, in the gravitino diagram 
with the quartic ipfi-f-cfyAi, vertex exists for the squark- 
gravitino production. 

Figure [5] shows squared matrix elements of the above 
two processes, qg — > qG (a) and gg — > gG (b), at the 
partonic CM energy VI — 2 TeV and at the partonic 
scattering angle cos^ = 0.5 as a function of the gravitino 



Diagrams by MadGraph u g -> ul gid 




graph 3 graph 4 

Fig. 1. Feynman diagrams for associated goldstino produc- 
tions with a squark, qg — >■ qG (top), and with a gluino, 
gg — > gG (bottom), generated by MadGraph. ul, go, and gld 
denote a sup, a gluino, and a goldstino, respectively. 

mass (7713/2). The squark and gluino masses are fixed at 1 
TeV. In the region of the small gravitino mass, < 
0.05, or in the high energy region ^13/2, the ratio 

of the squared matrix elements, |MgoidstinoP/|-MgravitinoP, 
is nearly unity, that is, both amplitudes agree well each 
other. We note that the squared matrix elements for the 
associated gravitino productions are proportional to "m'J^ 
as clearly seen in the log-log plot in Fig. [21 Associated 
productions of gravitino and gluino for the qq initial state 

q + q^~g + G, (8) 
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gravitino mass [GeV] 

Fig. 2. Squared matrix elements of gravitino (solid) and gold- 
stino (dashed) productions associated with a squark, qg — )■ qG 
(a), and with a gluino for the gg initial state, gg — >■ gG (b), 
at \fi = 2 TeV and cos^ = 0.5 as a function of the gravitino 
mass, where the squark and gluino masses are fixed at 1 TeV. 
The ratios of the squared matrix elements are also shown. 



which involves the ip-f-(f> and ip-X-A/^i vertices, can be also 
tested, as well as the above processes with an extra parton 
and crossed processes. 

Before turning to sample results, we also note that, in 
addition to the goldstino-gravitino equivalence test, the 
code was checked carefully by comparing with the ana- 
lytical squared matrix elements of eqs. (3), (5) and (7) 
in ref. [TT] for the above three partonic processes. The 
hadronic total cross sections at the SUSY benchmark points 
SPS7 and SPSS in figs. 8 and 11 in ref. [U can also be 
reproduced with the help of ME [5l[6]. The test by using 
the gauge invariance of the amplitudes is also mentioned 
in App. IA.4I 



3 Sample results 

In this section, we present some sample numerical results, 
using the new HELAS subroutines, which are presented in 
Appendix 1X1 and the modified MG, which is described in 
Appendix |B] 

In gauge-mediated SUSY breaking scenarios, the grav- 
itino is often the lightest supersymmetric particle (LSP), 
and its phenomenology depends on what the next-to-lightest 
supersymmetric particle (NLSP) is. While the lightest neu- 
tralino and the lighter stau are often the NLSP in minimal 
models of gauge mediation, gluinos can also be the NLSP, 
e.g., in split SUSY models and general gauge mediation 
models; see review papers pi ll2l and references therein. 

If gluinos are the NLSP and light enough, those pro- 
ductions can be explored in the early LHC data as well as 
in the Tevatron, and several studies have been performed 
for hadroproductions of a gravitino with a gluino (or a 
squark) [T^[H1[TT] . which lead to characteristic signals 
of monojet plus missing energy when a produced gluino 
(squark) promptly decays into a gluon (quark) and a LSP 
gravitino. We consider such scenarios for sample results of 
our code. 

Figure [3] presents total cross sections of each subpro- 
cess of associated gravitino productions with a gluino, 

pp/pp gG, (9) 

at the Tevatron-1.96TeV/LHC-7TeVfor the gravitino mass 
7713/2 = 10^^'' GeV as a function of the gluino mass. The 
masses of the left-handed and right-handed squarks, which 
appear in the t- and M-channel propagators, are fixed at 
1.5 TeV (dashed lines) and 2mg (dotted lines) for the qq 
subprocesses. The cross sections of associated productions 
with a squark, 

pp/pp and qG, (10) 

as a function of the squark mass are also shown in Fig. 2] 
for reference, where the productions of the left-handed 
and right-handed squarks are summed and their masses 
are taken to be same. The gluino mass is fixed at 1.5 
TeV (dashed lines) and 2mq (dotted lines). The CTEQ6L1 
parton distribution functions [TS] are employed, and the 
renormalization and factorization scales are fixed at the 
average mass of the final state particles, fXR = = 
(™g.i3 + '773/2)/2 mg,q/2. The cross section in the partial 
width r^(g)^q(g^G region is shown with a thin 

dotted line in Fig[31|l]), which will be discussed later. 

The major features of the production cross sections are 
following: 

— The cross sections of all the subprocesses scale with 
^^^2, that is, lighter gravitinos enhance the monojet 
signals, which can be interpreted as the direct lower 
bound for the gravitino massQ Note that, on the other 

^ The current lower bound from the Tevatron with the inte- 
grated luminosity of 87 pb~^ is about 10~^* GeV [TS], where 
all the SUSY particles except gravitino are assumed to be too 
heavy to be produced on-shell 117]. 
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Fig. 3. Total cross sections of each subprocess of associ- 
ated gravitino productions with a gluino, pp/pp — >■ gG, at 
the Tevatron-1.96TeV/LHC-7TeV for mg/a = 10"" GeV as 
a function of the gluino mass. The squark masses are fixed at 
1.5 TeV (dashed) and 2mg (dotted) for the qq subprocesses, 
where the cross section in the r,^ q > mq/2 region is shown 



with a thin dotted line. 



hand, the dijet signals which are produced through 
gluino-pair productions do not depend on the grav- 
itino mass. 

The gluino-gravitino associated productions through 
the gg initial state depend only on the gluino mass once 
the gravitino mass is fixed, while the 55-initiated pro- 
duction cross section depends not only on the gluino 
mass but also on the t- and u-channel-exchanged squark 
masses. It should be noted that those contributions 
are not decoupled in the large squark mass limit, and 
the heavier squark exchange increases the cross section 
since the goldstino-quark-squark couplings are propor- 
tional to the squark mass squared. Therefore, for heavy 
squark cases, the cross section of the qq channel can 
be larger than that of the gg channel even for the pp 
colhder, LHC. 
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Fig. 4. Total cross sections of associated gravitino productions 
with a squark, pp/pp — > qG, at the Tevatron-1.96TeV/LHC- 
7TeV for 77x3/2 = 10"^^ GeV as a function of the squark mass. 
The gluino masses are fixed at 1.5 TeV (dashed) and 2niq (dot- 
ted), where the cross section in the F-^^q > mg/2 region is 
shown with a thin dotted line. 



a gravitino is given by 



^9(9)- 



►9(9)G 



487rMpirn2 



(11) 



3/2 



Therefore, as the gravitino mass becomes small and/or 
the gluino (squark) mass becomes large, the partial width 
rapidly grows. For instance, in the 777.3/2 = 10""'^'^ GeV 
case, r = 1.12 x lO^^ GeV for TTigj,-) = 100 GeV, while 
r = 850 GeV for 777^(5-) =1.5 TeV. When the width be- 
comes a significant fraction of the mass, the goldstino cou- 
pling becomes strong and the perturbative calculations 
are unreliable. Thin dotted lines in Figs. |3] and |4] show 
such parameter regions. Moreover, the decay branching 
ratios depend on the entire SUSY mass spectrum; see, 
e.g., ref. [TT]. Therefore, the comprehensive analyses are 
necessary for the collider signatures such as monojet plus 
missing energy. 



Similar to the gluino productions through the qq sub- 
processes, the squark-gravitino associated productions, 
which are only produced through the qg subprocesses, 
involve the ^-channel gluino exchange diagram as well 
as the squark exchange one (see fig. H^top)), and the 
heavy gluino increases the cross section. 



Finally, it is worth noting that the partial width of a 
gluino (squark) decay into a gluon (massless quark) and 



4 Summary 

In this paper, we have added new HELAS subroutines to 
calculate helicity amplitudes with goldstinos, which ap- 
pear as the longitudinal modes of massive gravitinos in 
high energy processes, into the HELAS library. They are 
coded in such a way that arbitrary amplitudes with exter- 
nal goldstinos can be generated automatically by MadGraph, 
after slight modifications. All the codes have been tested 
carefully by making use of the goldstino equivalence the- 
orem and the gauge invariance of the helicity amplitudes. 
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Vertex 


Inputs 


Output 


Subroutine 


FFS 


FFS 


A in T~il i "i" 1 1 H p 


lORSGX, 


IROSGX 




FS 


F 


FSORGX, 


FSIRGX 




FF 


S 


HIORGX, 


HIROGX 


FFV 


FFV 


Amplitude 


lORVGX, 


IROVGX 




FV 


F 


FVORGX, 


FVIRGX 




FF 


V 


JIORGX, 


JIRQGX 


FFVV 


FFVV 


Amplitude 


IQRVVG, 


IROVVG 




FVV 


F 


FVVORG, 


FVVIRG 




FFV 


V 


J VI ORG, 


JVIROG 



Table 1. List of the new vertex subroutines in HELAS system. 



Total cross sections of associated gravitino productions 
with a gluino and a squark are also presented for the Teva- 
tron and the LHC as sample results. 
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A HELAS subroutines for goldstinos 

In this appendix, we list the contents of all the new HELAS 
subroutines that are needed to evaluate processes with an 
external goldstino based on the effective Lagrangian of 

In Apps. lA.ll to IA.3[ we explain vertex subroutines 
listed in Table [TJ which compute interactions among a 
goldstino with SM and SUSY particles. Those subrou- 
tines which we do not present in this paper are a FFSS 
type vertex, the last term in ((T]) . This contributes to some 
phenomenology in cosmology, e.g., the goldstino emission 
rate through qq — > gG and qg — qG [3], but it is mostly 
irrelevant in collider signatures. In App. IA.4I we briefly 
mention a test of our new subroutines by making use of 
the gauge invariance. 

Before turning to vertex subroutines, we note that 
the existing subroutines for flowing-in (IXXXXX) and -out 
(OXXXXX) spinor wavefunctions can be applied for a gold- 
stino without any modification since it is just a spin-1/2 
particle. 



A.l FFS vertex 

The FFS vertices involving a goldstino are obtained from 
the interaction Lagrangian among a fermion (/), a gold- 



Vertex 


New subroutine 


Reference subroutine 


FFS 


IDRSGX, IROSGX 


lOSXXX 




FSORGX 


FSOXXX 




FSIRGX 


FSIXXX 




HIORGX, HIROGX 


HIOXXX 



Table 2. Reference HELAS subroutines in the present HELAS 
library 3 to the new ones. 



stino and a scalar boson (S): 

rpFs = ^(GC(I)Pl + GC(2)Ffl)/ S* + h.c. 



(12) 



with the chiral-projection operator Pj^/l = ■^(1 i 75). 
The interaction vertices are exactly same as those in the 
present HELAS library, i.e., the Yukawa- type vertex. There- 
fore, we simply reuse the existing lOSXXX, FSOXXX and 
HIOXXX subroutines for lORSGX, FSORGX and HIORGX, re- 
spectively, while IROSGX, FSIRGX and HIROGX need cou- 
pling modification due to the Hermitian conjugate as 



GC'(l) = (GC(2))*, 
GC'(2) = (GC(1))*, 



(13) 



in order to reuse lOSXXX, FSIXXX and HIOXXX, respectively. 
The correspondence between the new HELAS subroutines 
and the existing ones are shown in Table [51 see also the 
HELAS manual [3" for details of these subroutines. 

GC ( 1 ) and GC (2) are the relevant left- and right-coupling 
constants. For instance, in the case of the quark-goldstino- 
squark interaction without squark mixing, q-G-qa, those 
couplings read 




= GFFSL(l) = GFFS 
= GFFSL(2) = 



GC(1) 
GC(2) 



GFFSR(l) 
GFFSR(2) 







-GFFS 



for a = L, (14) 
for a = R, (15) 



with 



GFFS = im? /V3M 



PI W3/2 



(16) 



in the quark massless limit. We note that, in the HELAS 
convention, the factors of i in the goldstino coupling con- 
stants are necessary to give a correct phase between the 
amplitudes which involve particles in minimal supersym- 
metric standard model (MSSM) plus goldstinos. 



A. 2 FFV vertex 

The FFV vertices involving a goldstino are obtained from 
the interaction Lagrangian among a fermion, a goldstino 
fermion and a vector boson [V^): 



Cff, = -iVi[7^7l(GC(l)PL +GC(2)P^)/9^F; 



h.c. 

(17) 



Although both a goldstino and a gaugino are Majorana 
particles in most cases, the Hermitian conjugate term is 



6 



K. Mawatari and Y. Takaesu: HELAS and MadGraph with goldstinos 



necessary for MG; practically, either the first or second term 
is used in calculations of amplitudes with the coupling 
constant 



GC(1) = GC(2) = GFFV 



where the inputs FMASS and FWIDTH are the mass and the 
width of the ofF-shell fermion, mp and Fp. The output 
FVORG (6) is a complex six-dimensional array and gives the 
ofF-shell fermion wavefunction multiplied by the fermion 
i//2\/6-/Vfpi TO3/2 (18) propagator and its momentum as 



in units of GeV ^, where m/ is the gaugino mass. 

We note that the input and output parameters of the 
subroutines for the FFV goldstino vertices are completely 
same as those of the existing FFV-type subroutines, though 
the vertex structure computed in the subroutines is dif- 
ferent. 



A. 2.1 IDRVGX 

This subroutine computes an amplitude of the FFV gold- 
stino vertex from wavefunctions of a flowing-in fermion, a 
flowing-out goldstino and a vector boson, and should be 
called as 

CALL IDRVGX (FI,FD,VC,GC , VERTEX). 

The inputs FI(6) and FD(6) are complex six-dimensional 
arrays which consist of the wavefunction and the four- 
momentum of the flowing-In and -Out Fermion as 

= (SReFI(5),5ReFI(6),3mFl(6),5TOFl(5)), 
p'^ = (3fieF0(5),5ReFD(6),3mF0(6),3TOFD(5)), 

respectively. VC (6) is a complex six-dimensional array which 
contains the Vector boson wavefunction and its momen- 
tum as 

= (3fieVC(5),3fieVC(6),5TOVC(6),3mVC(5)). 

The input GC(2) is the coupling constant in ([T51) . The 
output VERTEX is a complex number: 

VERTEX = (FD)[^, f ](GC(1)Pl + GC(2)Ffl)(Fl), (19) 

where we use the notations 



(FVORG) = (F0)[^, f ](iGC(1)Pl + iGC{2)Pr) 

+ mp) 



— rrip + impFp 



and 



FV0RG(5) = F0(5) + VC(5), 
FV0RG(6) = F0(5) + VC(6), 

where we use the notation 



(23) 



(24) 
(25) 



(FVORG) =(FV0RG(1), FV0RG(2), FV0RG(3), FV0RG(4)), 

(26) 

and the momentum p is 

p^ = (5ReFVORG(5), 5ReFVORG(6), 

3toFV0RG(6), 3toFV0RG(5)). 

A.2.3 JIORGX 

This subroutine computes an off-shell vector current made 
from the interaction of a flowing-in fermion and a flowing- 
out goldstino by the FFV goldstino vertex, and should be 
called as 

CALL JIORGX (FI,FO,GC,VMASS,VWIDTH , JIORG), 

where the inputs VMASS and VWIDTH are the mass and 
the width of the vector boson, my and Fy- The output 
JIORG (6) is a complex six-dimensional array and gives 
the off-shell vector current multiplied by the vector boson 
propagator and its momentum as 



JIORG(/^+ 1) 



-9 



(FI) = 



FI(1)N 
_ I FI(2) 



FI(3) ' 



vFl(4)y 

(FO) = (F0(1),F0(2),F0(3),F0(4)), 
1/^ =VC(/i + l). 



(20) 

(21) 
(22) 



X (F0)[^,7,](iGC(l)Pi + iGC(2)PH)(Fl) (27) 
for the massive vector boson, or 

JIORG(Ai+ 1) = ^ 
q^ 

X (F0)[^,7'^](zGC(l)Pi + iGC(2)Pfl)(Fl) (28) 
for the massless vector boson, and 



A. 2. 2 FVORGX 

This subroutine computes an off-shell fermion wavefunc- 
tion made from the interaction of a vector boson and a 
flowing-out goldstino by the FFV goldstino vertex, and 
should be called as 

CALL FVORGX (FO,VC,GC, FMASS, FWIDTH , FVORG), 



JI0RG(5) = -FI(5) + F0(5), 
JI0RG(6) = -FI(6) + F0(6). 



(29) 
(30) 



Here, q is the momentum of the off-shell vector boson, 

=(5ReJIORG(5), 3?eJIORG(6), 5mJI0RG(6), 3mJI0RG(5)). 

Note that we use the unitary gauge for the massive vector 
boson propagator and the Feynman gauge for the massless 
one, according to the HELAS convention [3]. 
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A. 2. 4 IRDVGX 

This is essentially the same subroutine as IDRVGX but 
modifies the input coupling GC as in (jl3p since this com- 
putes the Hermitian conjugate vertex of lORVGX. 



A. 2. 5 FVIRGX 

This subroutine computes an off-shell fermion wavefunc- 
tion made from the interaction of a vector boson and 
a flowing-in goldstino by the FFV goldstino vertex, and 
should be called as 

CALL FVIRGX (FI,VC,GC,FMASS,FWIDTH , FVIRG). 

What we compute here is 



(FVIRG). -J£±i^£l 



— TTip + imprp 
X (iGC(l)*P« + iQC{2)*PLM f ](FI), (31) 



and 



(32) 
(33) 



(34) 



FVIRG(5) = FI(5) - VC(5), 
FVIRG(6) FI(6) - VC(6), 

where we use the notation 

(FVIRG(1)\ 
FVIRG(4)/ 

and the momentum p is 

= (5ReFVIRG(5),5ReFVIRG(6), 
3mFVIRG(6), 3mFVIRG(5)). 



A.2.6 JIROGX 

This is essentially the same subroutine as JIDRGX but 
modifies the input coupling GC as in (jl3p since this com- 
putes the Hermitian conjugate vertex of JIORGX. 



coupling constant and the gauge coupling constant of the 
involving gauge boson. For instance, in the case of the 
goldstino-gluino-gluon-gluon interaction, G-g-g-g, those cou- 
plings are 

GC(1) = GC(2) = GFFV * Gs = -imfg,/2V6Jlpi 1713/2 

(36) 

in units of GeV~^, where Gs = g^ is the strong coupling 
constant. 



A.3.1 lORVVG 

This subroutine computes an amplitude of the FFVV gold- 
stino vertex from a flowing-in fermion, a flowing-out gold- 
stino and two vector bosons, and should be called as 

CALL IORVVG(FI,FO,VA,VB,GC , VERTEX). 

What we compute here is 

VERTEX = (FO)[f ^ f ''](GC(1)Pl + GC(2)P7j)(Fl), (37) 
where we use the notations 



= Mk{n + 1) 

v''''' = min + 1) 



(38) 
(39) 



A.3.2 FVVORG 

This subroutine computes an ofF-shcll fermion wavefunc- 
tion made from the interaction of two vector bosons and 
a flowing-out goldstino by the FFVV goldstino vertex, and 
should be called as 

CALL FVVORG (FD,VA,VB,GC,FMASS,FWIDTH , FVVORG). 
What we compute here is 

(FVVORG) = (FO)[f f ^](iGC(l)PL + iGC(2)Pfl,) 

p"^ — m|n -I- imprp ' 



A.3 FFVV vertex 



and 



The FFVV vertices involving a goldstino are obtained from 
the interaction Lagrangian among a fermion, a goldstino 
and two vector bosons: 

-Cffw = ^r''^v5[7^7"](Gc(l)PL + Gc(2)P;^)r 

-l-h.c. (35) 

with the structure constant Z"'"^, which can be handled 
by the MG automatically. As in the FFV vertex case, we 
need the Hermitian conjugate term although a goldstino 
and a gaugino are Majorana particles in most cases. The 
coupling constant GC is the product of the FFV goldstino 



FVV0RG(5) = F0(5) + VA(5) + VB(5), (41) 
FVV0RG(6) = F0(6) + VA(6) + VB(6), (42) 

where we use the notation 

(FVVORG) =(FVV0RG(1), FVV0RG(2), FVV0RG(3), FVV0RG(4)), 

(43) 

and the momentum p is 

= (5ReFVVORG(5), 3?eFVVORG(6), 
3mFVV0RG(6), 3mFVV0RG(5)). 
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A. 3. 3 JVIDRG 



where we use the notation 



This subroutine computes an off-shell vector current made 
from the interaction of a vector boson, a flowing-in fermion 
and a flowing-out goldstino by the FFVV goldstino vertex, 
and should be called as 

CALL JVIORG(FI,FO,VC,GC,VMASS,VWIDTH , JVIORG). 

What we compute here is 



JVIORG(^+ 1) = 



q- — niy + irnvTv 



m 



V 



X (F0)[7„ f ](zGC(1)Pl + zGC(2)Pi^)(Fl) (44) 
for the massive vector boson, or 



JVIORG(/i+ 1) = — 

q 

X (FD)[7^f](iGC(l)Pi- 
for the massless vector boson, and 



iGC(2)PR)(Fl) (45) 



JVI0RG(5) = -FI(5) + F0(5) + VC(5), 
JVI0RG(6) = -FI(6) -I- FD(6) -f VC(6), 



(46) 
(47) 



where momentum q is 



(3fieJVIORG(5), 3?eJVIDRG(6), 
SmJVI0RG(6), 3mJVI0RG(5)). 



A.3.4 IROVVG 



This is essentially the same subroutine as lORVVG but 
modifies the input coupling GC as in psp since this com- 
putes the Hermitian conjugate vertex of lORVVG. 



A. 3. 5 FVVIRG 

This subroutine computes an off-shell fermion wavefunc- 
tion made from the interaction of two vector bosons and 
a flowing-in goldstino by the FFVV goldstino vertex, and 
should be called as 

CALL FVVIRG (FI,VA,VB,GC,FMASS,FWIDTH , FVVIRG). 

What we compute here is 



(FVVIRG = — ^A, ^ 



and 



— rrip + impFp 
X {iGC{l)*PR + iGC{2)* PlW\ f^]{7l), (48) 



FVVIRG(5) = FI(5) - VA(5) - VB(5), (49) 
FVVIRG(6) = FI(6) - VA(6) - VB(6), (50) 



(FVVIRG) = 



and the momentum p is 



/FVVIRG(1)\ 
FVVIRG(2) I 
FVVIRG(3) ' 

\FVVIRG(4)/ 



(51) 



= (5ReFVVIRG(5), 3fieFVVIRG(6), 
5mFVVIRG(6), 3mFVVIRG(5)). 



A.3.6 JVIROG 

This is essentially the same subroutine as JVIORG but 
modifies the input coupling GC as in p3p since this com- 
putes the Hermitian conjugate vertex of JVIORG. 



A. 4 Checking for the new HELAS subroutines 

In addition to the goldstino-gravitino equivalence test in 
Sect. [21 the new HELAS subroutines are tested by using the 
gauge invariance. Helicity amplitudes involving external 
gluons are expressed as 



Mx^=e^{pg,Xg)T^^ 



(52) 



by extracting one of the gluon wavefunctions, £^{pg, Xg), 
with a momentum pg and a helicity Xg. The identity for 
the SU{3) gauge invariance 



(53) 



can be used for checking HELAS subroutines. In particular, 
we test the invariance of the following processes; 



qg qG 

gg gGg 



for FFS and FFV subroutines, 
for FFV and FFVV subroutines. 



We also test the agreement of the helicity-summed ampli- 
tude squared evaluated in arbitrary Lorentz frames. 



B Implementation of goldstinos into 
MadGraph 

In this appendix, we describe how we implement goldsti- 
nos and their interactions into MG. 

First, using the default mssm model in MG/MEv4 0, we 
make our new model directory, mssm_goldstino, where 
we define a goldstino (particles.dat) and its interac- 
tions with SM and SUSY particles (interactions.dat 
and couplings . f); we show the coupling constants for 
each goldstino vertex involving SUSY QCD particles in 
Table [3] as examples^ Then we add all the new HELAS 



^ It should be noted that the effective Lagrangian ([T| allow 
to calculate the FFV and FFVV interactions only for massless 
gauge bosons. 
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3-point couplings GC 



FFS 


q gld 


qi 


GFFSL 




q gld 




GFFSR 


FFV 


go gld 


g 


GFFV 


4-point 


couplings 




GC 


FFVV 


go gld 


g g 


GFFVV = GFFV*Gs 



Table 3. List of the coupling constants for each goldstino ver- 
tex involving SUSY QCD particles. All the particles and the 
coupling constants are written in the MG notation, gld stands 
for a goldstino, q represents a light quark, and ql/qr is a 
left/right-handed squark. g and go are a gluon and a gluino, re- 
spectively. GC is a coupling constant defined in each subroutine 
in Add. 1X1 



15. J. Pumplin, D. R. Stump, J. Huston, H. L. Lai, 
P. M. Nadolsky and W. K. Tung, JHEP 0207 (2002) 012. 

16. A. A. Affolder et al. [CDF Collaboration], Phys. Rev. Lett. 
85 (2000) 1378. 

17. A. Brignole, F. Feruglio, M. L. Mangano and F. Zwirner, 
Nucl. Phys. B 526 (1998) 136 [Erratum-ibid. B 582 (2000) 
759]. 



subroutines for goldstinos to the HELAS library in MG. Al- 
though the present MG can handle spin-1/2 particles, we 
simply replace spin-3/2 gravitinos in the recent modified 
MG [2| by goldstinos since the FFV goldstino vertex is non- 
renormalizable and has to be distinguished by the exist- 
ing renormalizable FFV vertex and the flowing-in and -out 
goldstinos require different subroutines as in the gravitino 
case. We note that, in order to achieve the above simple 
modification, we put "G" into the end of names of each 
gravitino subroutine like 

lORSXX IDRSGX. 
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